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A NEW UPPER BOUND FOR FINITE ADDITIVE BASES 


C. SiNAN GUNTURK AND MELVYN B. NATHANSON 

Abstract. Let n(2, k) denote the largest integer n for which there exists a set 
A of k nonnegative integers such that the sumset 2A contains {0,1, 2,..., n — 
1}. A classical problem in additive number theory is to find an upper bound 
for n(2,k). In this paper it is proved that limsupj,_,j^ n{2,k)/k‘^ < 0.4789. 


1. An EXTREMAL PROBLEM FOR FINITE BASES 

Let No and Z denote the nonnegative integers and integers, respectively, and let 
|A| denote the cardinality of the set A. 

Let A be a set of integers, and consider the sumset 

2A = {a-\- a : a,a € A}. 

Let ^ be a set of integers. The set A is a basis of order 2 for S' if S' C 2A. 
The set A is called a basis of order 2 for n if the sumset 2A contains the first n 
nonnegative integers, that is, if A is a basis of order 2 for the interval of integers 
[0, n — 1] := {0,1, 2,..., n — 1}. We define n(2. A) as the largest integer n such that 
A is a basis of order 2 for n, that is, 

n(2. A) = max{n : [0, n — 1] C 2A}. 

Rohrbach introduced the extremal problem of determining the largest integer 
n for which there exists a set A consisting of at most k nonnegative integers such 
that A is a basis of order 2 for n. Let 

n(2, k) = max{n(2. A) : A C No and | A| = k}. 

Rohrbach’s problem is to compute or estimate the extremal function n(2,k). The 
set A is called an extremal k-basis of order 2 if |A| < k and n(2. A) = n(2, k). 

For example, n(2,1) = 1 and n(2, 2) = 3. The unique extremal 1-basis of order 
2 is {0}, and the unique extremal 2-basis of order 2 is {0,1}. For fe = 3 we have 
n(2, 3) = 5, and the extremal 3-bases of order 2 are {0,1, 2} and {0,1, 3}. If fc > 2 
and A is an extremal fc-basis of order 2, then 0,1 S A. If A is a finite set of k 
nonnegative integers and n(2. A) = n, then n ^ A. If a G A and a > n, then the set 
A' = (A \ {a}) U {n} has cardinality k, and n(2. A') > n -|- 1 > n(2. A). Therefore, 
if A is an extremal fc-basis of order 2 and n(2, fe) = n, then 

{0,1} C A C {0,l,2,...,n- 1} C 2A. 
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If A is an extremal fc-basis for n, then |A| = k and A C {0,1, 2,..., n — 1}. 

Rohrbach determined order of magnitude of n(2, k). He observed that if H is a 
set of cardinality k, then there are exactly ordered pairs of the form {a, a') 

with a, a' & A and a < a'. This gives the upper bound 

n{2,k)<(^^^^'^=^ + 0{k). 

To derive a lower bound, he set r = [k/2] and constructed the set 
H = {0,1, 2,..., r - 1, r, 2r, 3r,..., (r - l)r}. 

We have 

|H| = 2r - 1 < fc 


and { 0 , 1 ,.. 

., r^} C 2A. Then 


n( 2 ,H)>r^ + l> ^ +l = - 

and so 

n{2,k) > ^ + 0{k). 

Thus, 

li„i„f"P;'')> 1 = 0.25 

n^oo 4 

and 

n(2,k) 1 

limsup ——— < 7 : = 0.5. 

n— *00 k 2 


It is a open problem to compute these upper and lower limits. Mrose urn proved 
that 

liminf > ^ = 0.2857..., 

n—»oo 7 

and this is still the best lower bound. Rohrbach used a combinatorial argument to 
get the nontrivial upper bound 

limsup^^^^l^ < 0.4992. 

Moser |3] introduced Fourier series to obtain 

limsup < 0.4903, 

and subsequent improvements by Moser, Pounder, and Riddell 0] produced 

limsup < 0.4847. 

Combining Moser’s analytic method and Rohrbach’s combinatorial technique, Klotz | 2 ] 
proved that 

limsup "^^’^^ < 0.4802. 

In this paper, we use Fourier series for functions of two variables to obtain 

limsup< 0.4789. 

n—*oo rZ 
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We note that Rohrbach used a slightly different function n{2,k): He defined 
n(2, k) as the largest integer n for which there exists a set A consisting of fc + 1 
nonnegative integers such that the sumset 2A contains the first n + 1 nonnegative 
integers. Of course, Rohrbach’s function and our function have the same asymp¬ 
totics. 


2. Moser’s application of Fourier series 


In this section we describe Moser’s use of harmonic analysis to obtain an upper 
bound for n(2, k). Let A be an extremal fc-basis of order 2. Let r 2 ,A{j) denote the 
number of representations of j as a sum of two elements of A, that is, 

i'2,a 0’) = card ({(oi, 02 ) € A x A ai + 02 = j and oi < 02 }). 

We introduce the generating function 

aGA 

k = fA{l) = \A\ 

fA{qf + fA{q^) 


Then 

and 


= X] '^2,AU)q^- 

jG2A 

If [0, n — 1] C 2A, then r 2 ,A{j) > 1 for all 0 < j < n — I. Hence there exist integers 
S{j) > 0 such that 

fA{q? + fA{q^) 


where 


Let 


Sij) = 


= l + q + q^ + ... + q--^+ W, 

jG2A 

C 2 .a 0')-1 if j e {0,l,...,n-1}, 
C 2 ,a 0 ) otherwise. 


A(g) = ^ 5{j)q^. 

jG2A 


Then A{q) > 0 for q > 0, and 

/A(g)^ + /A(g^) 


= l + q + q^ + --- + q 


n — 1 


A{q). 


Evaluating the generating function identity o at g = 1, we obtain 

k'^ + k 


( 2 ) 

Since A(I) > 0, we have 


= n + A(l). 


n < y -HO(fc). 


The strategy is to find a lower bound for A(l) of the form 

A(l) > ck^ + 0{k) 
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for some c > 0, and deduce 

n < ^ 0{k). 

We obtain a simple combinatorial lower bound for A(l) by noting that if Oi, 02 G 
A and nj2 < ai < a 2 , then 01+02 > n. Let i denote the number of elements a € A 
such that o > n/2. Then 


( 3 ) 

Let 


^(1) > = XI’"2^0) > —5— > y- 


j>n 


j>n 


u = 


be a primitive nth root of unity. Let r be an integer not divisible by n. Then 

1 + w’’+ + ... + = 0 

and so 

/a(w'’)^ + /a(w^’') _ 1 , , .r , , .2r , , , .(n-l)r 


= 1 + W’’ + + • • • + W 




Applying the triangle inequality, we obtain 

I/aK)"+/a(c.2o| ^ \fA{con\^-k 


A(l) > |A(u;’-)| = 


> 


Let 

Then 

( 4 ) 

and 

( 5 ) 


2 - 2 
M = max{|/^(w'’)| : r ^ 0 (mod n)}. 

0< M <k 


We can also obtain an analytic lower bound for A(l). For all integers r not 
divisible by n, we have 




^ ^ ^2-Kirajn 

= 

XI cos(27rro/n) + i sin(27rra/n) 

aGA 


aeA 


and so 


cos(27rro/r 


aeA 


and 


< M 


< M. 


XI sin(27rro/?' 
a^A 

Let ip(t) be a function with period 1 and with a Fourier series 

00 00 

ip{t) = r cos(27rrt) + sin(27rrt) 


r=0 
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whose Fourier coefficients converge absolutely, that is, 


1“'’! + X! \br\ < OO. 
r—0 r—1 


Let 




r = 0 
n|r 


For any integer a we have 

OO OO 

= EE®’- cos( 27 rra/n) + EE br sin( 27 rra/r 


aGA 


a£A r—0 


a^A r—0 


= E“’-E cos( 27 rra/n) + E^’’ E sin( 27 rra/n) 

r=0 a£A r—1 a^A 

OO OO OO 

= E“-E cos( 27 rra/n) + sin( 27 rra/n) + A: , 

>■=0 aeA 


'■=1 aeA 


n|r 


and so 




aeA 


<M^(|a,| + |6,|) + fcC. 


Let ai and 02 be real numbers such that 

ip(t) > ai for 0 < t < 1/2 

and 

(p{t) > 02 for 1/2 < t < 1. 

Recall that i denotes the number of elements a G A such that n /2 < a < n — 1 . 
Then 

ip ^ > {k — t)a\ + la 2 = kai — (ai — 02)^. 

aSA 

We obtain the inequality 


( 6 ) 


kai — (ai — 02)^ < M y](|ar-| + |A'r|) + kC. 


r=0 

n-fr 


In this way, the function p{t) produces a lower bound for M, which, by @, gives 
a lower bound for A(l). 

Moser applied inequality to the function 

p{t) = t cos(47rt) + sin( 27 rt), 

whose nonzero Fourier coefficients are 02 = 1/2 and 61 = 1 . Then C = 0 for n > 3 , 
and 

3 M 

I < 

^n/ 




aeA 
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The function ip(t) satisfies the inequality 


^ J 5 for 0 < t < 1/2 


and so 


This implies that 


EH/) 


> 


k-e u k-u 


aeA 


M>- 
- 3 


EH/) 


aeA 


> 


k — 4:i 


and we obtain the analytic lower bound 

A(l) > fcill! - 

^ ^ - 18 2 

Recalling the combinatorial lower bound © 

M'^-k 


we obtain 


A(l) > max 


A(l) > 

{k-Uf k fc2 k 

18 1 2 “ 2' 


Inserting this into inequality we obtain 

+ fc 


, , k"^ k 

= n + A(l)>nH-, 

2 ^ - 98 2 ’ 


and so 


n < 


1 1 
2 ~ ^ 


+ fc < 0.4898fc^ + k. 


3. Fourier series in two variables 

We shall modify Moser’s method to obtain a better lower bound for A(l). We 
use the same notation as in the previous section. In particular, i denotes the 
number of integers a G A such that a > n/2. Let L denote the number of pairs 
(oi, 02 ) G A X A such that Oi + 02 > n. Then L > and k^ — L is the number of 
pairs (oi, 02 ) G Ax A such that oi + 02 < n — 1. We have the combinatorial lower 
bound 

(7) A(l)>^r2.AH = ^>f 

j>n 

Let ip{ti,t 2 ) be a function with period 1 in each variable and with a Fourier 
series 

whose Fourier coefficients converge absolutely, that is, 

EE |‘^(l'l,'I’ 2 )| < 00. 

ri GZ r2€Z 
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We choose with zero mean, that is, 

(, 5 ( 0 , 0 )= f f (f{ti,t 2 )dtidt 2 = 0 . 

Jo Jo 

Let 

-Ri = {(ill ^ 2 ) G [0i 1) X [0,1) : ti + ^2 < 1} 

and let 

i?2 = {(ill ^ 2 ) G [0i 1) X [0,1) : ii + i2 > 1} 

If 01,02 G A and 01+02 < n — 1 , then (oi/n, 02/n) G Ri. If oi + 02 > n, then 
{aijn, 02/n) € i?2- 

Let «! and a2 be real numbers such that 

+(iiii 2 )>ai for (ti,t 2 ) G i?i 

and 

+(ii,i 2 )>a 2 for (ti,t 2 ) G R 2 - 

We choose the function ip{ti,t2) such that 

cxi > q;2 - 

Then 

( 8 ) if f —, —) > (A:^ — L)ai + La2 = aik‘^ — {ai — a2)L. 

^ , \ n n / 

CLi CL2^A 

We can rewrite this sum as follows: 

E E ^(t’t) = E E e 

ai^Aa2GA aiGAa2^Ari^Zr2^'^ 

r\&r2& ai^A a2^A 

= X! X! ‘^(^1>^2 )/a(w’'^)/a(w’'"). 

Consider the partition of the integer lattice 1 ? = Rq U Ri U 5'2: 

Ro ={(7'1 iT2 ) G : Ti = 7-2 = 0 (mod n)} 

Si ={(ri,r2) G Z^ : ri = 0 (mod n),r2 ^ 0 (mod n)} 

U {(ti, r2) G Z^ : ri ^ 0 (mod n), r2 = 0 (mod n)} 

5'2 ={(Ti,r2) G Z^ : ri ^ 0 (mod n),r2 ^ 0 (mod n)}. 

We define Co, Ci, and C2 by 

Ci= ^ |(/>(ri,r2)|. 

(i’l.r2)eSi 

Recall that |/a (w'') | < M if r is not divisible by n and |/a (w’’) | < fc if r is divisible 
by n. Then 

EE.(=t) 

qi^A 0,2^a 


(9) 


< Cofc2 + CifcM + C2M2. 
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Combining inequalities 0 and we obtain 

aik^ — (ai — a2)L < Cok^ + CikM + C2M^. 
Since ai > a 2 , we have 

{ai - Co)k^ -CikM- 


L > 


We define 


0.1 — 02 
M 


M = 


Since 0 < M < k, we have 

0 < Ai < 1. 

By inequality 0 , we have 2A(1) > L, and so 


2A(1) ^ ^ (ai — Cq) — Cin — C2/i^ 


k'^ k'^ a\ — 02 

By inequality we also have 2A(1) > — k, and so 


( 10 ) 


2A(1) 

fc 2 


> max ( Ai , 


(ai - Co) - CiAi - C2Ai' 


1 


Oi — 02 

Since the series of Fourier coefficients of ip{t\,t 2 ) converges absolutely and since 
(^(0, 0) = 0, we can arrange the Fourier series in the form of a sum over concentric 
squares 

OO 

R=l max(|ri|,|r 2 |)=-R 

For any e > 0 there exists an integer N = N{e) such that 

OO 

E E |<^(i'i,i’2)| < e{oi - 02). 

n—N max(|ri \ ,\r 2 \)—n 

For all n > N, we shall approximate the sums Co,Ci, and C 2 by 0, Caxiai, and 
Cmain, respectively, where 


and 


Then 


Caxiai = Y + I<^(a0)|) 

rez 

Cmain = E E I <^(^1 >’’z) | ■ 


riez r 2 ez 

r2^0 


|(ai — Co) — CiAi — C2^J-^ — [oi — CaxialAi ~ CmainA*^) | 

= |Co + (Cl — Caxial)/i + (C 2 — Cinain)Ai^| 

< |Co| + |Ci - Caxiall + IC2 - CmainI 

< E \'f>{ri,r2)\ 

max(|ri|,|r2|)>Af 

< e(ai - q; 2 ), 
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and so 


(ai - Co) - CiAi - C2 m' 


dl — Q!2 

It follows from inequality m that 


^1 CaxialM CaaainM 

di — d2 


< e. 


2A(1) 


> max I /r , 


0^1 CaxialM CmainM 

dl — d 2 


^ k' 


Let 

( 11 ) 


. r / 2 ^axial/^ C^mainM 

p = mi max I p , — 


0<M<1 


ai — 0.2 


From TO and the definition of p in we now have 


2A(1) 1 

-p— >p-e--. 


Applying identity 0, we obtain 

+ fc 


= n + A(l) > n + 


{p — e)fc^ — k 


Therefore, 


n < 


1- p + e 


where the number p depends only on the function and e > 0 can be 

arbitrary small. 

It is clear that we always have p > 0, and that p > 0 if and only if di > 0. It 
is also clear that when di > 0, we have p = where ^ is the unique solution in 
[0,1] to the quadratic equation 

^2 _ Caxial^ Cmain^ 


O-l — 0-2 


i.e., 


(dl d 2 T Ciaiain)'^ “f CaxialC 0^1 — O 5 
which yields the formula 

“Caxial + \/ C'axial + 4di(di — d 2 + Cmain) 


( 12 ) 


P = 


2 dl(dl d 2 “t“ Cmain) 


Hence we have an optimization problem in which we maximize p over all real 
valued functions p defined on the unit square [0,1)^ such that p has zero mean and 
tp > 0 on i?i. We do not know the optimal function for this problem, but we have 
found a simple piecewise polynomial function that improves Klotz’s upper bound 
for n(2, k). Before we proceed to the main result of this paper, which also includes 
the definition of this function, let us present some of the heuristics which have lead 
us to our “educated guess.” 

First, without loss of generality, we may assume that di = 1. Note that we then 
necessarily have 

idi + id2< JJ ip{ti,t2)dtidt2 + jj ip{ti,t2)dtidt2 =0 
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SO that 02 < —1- We also have 


Caxial > 




{(p{0,t) + ip{t,0)) dt 


> 2, 


and 


C'rnain ^ 




f if{t, 1 — t)dt 

Jo 


> 1 . 


In any case we are interested in the positive root of the equation 

— 1 = 0 


where k = (1 — a 2 + Cmain) > 3 and r = Caxiai > 2. Clearly, the smaller k and r 
are, the larger this root will be. The bounds k > 3 and r > 2 already imply that 
^K,T < hence p = In reality, 02 < ~1 because equality can happen only 

if ip is constant on both Ri and i? 2 , in which case ip is not absolutely summable. 
This results in the heuristic that if we try to push 02 close to —1, then Caxiai and 
Cmain will become large, and conversely if we try to push Caxiai and Cmain close to 
their respective minimum values, then p may not be bounded from below on i ?2 
by a small value. The right trade-off between these two competing quantities will 
result in the solution of this optimization problem. 

It is interesting to note that the value of p is fairly robust with respect to varia¬ 
tions in K and r, which we will only be able to estimate but not compute exactly. 
The following lemma gives an explicit estimate for this purpose: 


Lemma 1. Let o,nd ^ko,to respective positive roots of the equations 

- 1 = 0 and -b - 1 = 0. Let p = and po = C^o.ro- 
min(K, kq) ^ 3 and min(T, tq) > 2, then 

(13) |p-po| < ^|K-Ko|-b^|r-To|. 

The proof of this lemma is given in the Appendix. Now we can state and prove 
the main theorem of this paper. 


Theorem 1. 


limsup "*'^^^^ < 0.4789. 

n —>00 


Proof. We define the function p{ti,t 2 ) on the unit square [0,1)^ by 

/i/,x f, i \ _ / li (C,^2)€i?i 

(14) <b>(ii,t2)-| (C,t2)ei?2. 

Then 

ai = 1 . 

Computation of the three other parameters used in formula OD for p yields 

«2 = 1 -^ = - 3 . 72470 ..., 


and 


2.90278 < Caxiai < 2.90289, 
4.75145 < Cmain < 4.76146. 
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Taking kq = 1 + 3.72471 + 4.76146 = 9.48617, and tq = 2.90289, we obtain po > 
0.04240, |k — Kol < 0.01002 and |t — tqI < 0.00011, so that |p—pol < 0.0002. Hence 

P > Po - Ip - Pol > 0.0422, 

and consequently choosing e sufficiently small, 

n < 0.4789 + 

The details of the computations are in the Appendix to this paper. This completes 
the proof. □ 


4. Open problems 


A major open problem concerning the extremal function 


n(2, k) = max{n(2. A) : A C No and | A| < k}. 


is to compute liminf„^oo 11 ( 2 , and limsup„^oo n(2, and to determine 

if the limit 


lim 

n^oo 


n(2, k) 


exists. We have no conjecture about the existence of this limit, nor about the values 
of the lim inf and lim sup . 

It is also difficult to compute the exact values of the function n(2, k). 
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We can generalize the extremal functions n{2,A) and n(2, k) as follows. Let A 
be a finite set of integers, and let m{2,A) denote the largest integer n such that 
the sumset 2A contains n consecutive integers. Let 

m(2, k) = max{m(2, A) : A C X and | < k}. 

Let i{2, A) denote the largest integer n such that the sumset 2A contains an arith¬ 
metic progression of length n, and let 

£(2, k) = max{.^(2, A) : A CZ and | A| < k}. 

We can also define the extremal function 


n'(2,k) = max{n(2, A) : A C Z and |A| < k}. 


Then 


n(2, A) < n'(2, A) < m(2, A) < 1(2, A) 


and so 


n(2, k) < n'{2, k) < m{2, k) < t{2, k). 

For any integer t and set A, we have the translation A + t = {a+ t -.aG A}. 
The functions i and m are translation invariant, that is, 1(2, A + t) = i(2. A) and 
m{2,A + t) = m(2,A). We also have the trivial upper bound £(2,k) < (^ 2 ^), but 
it is an open problem to obtain nontrivial upper bounds for any of the extremal 
functions n'{2, k), m{2, k), or £{2, k). 

Appendix 

We describe here the computations. 

Proof of Lemma m We start with the formula 


—T + + 4k 2 


2k t + + 4k 

We next evaluate the partial derivatives of with respect to k and r: 


(15) 





4 


(16) 


•s/r^ -b 4 k (t -|- ^/t'^ + 4 k)^ ' 
2 


•\/t^ -b 4 k (r -I- -s/r^ -b 4 k) 


from which it follows that in the set {(k,t) : k > 3, t > 2}, we have ^ 



hence the result of the lemma. 


□ 
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Absolute summability of ip. While the result explained in this subsection is 
elementary, we will provide a certain amount of detail in its derivation because our 
main concern is more than absolute summability of ip. We would like to provide 
explicit estimates on the rate of the convergence; this will be necessary later in 
the section when we will analyze the accuracy of the numerical computation of the 
constants domain and Coaxial ■ 


Lemma 2. Let f be a smooth function on [0,1]. Then for all L > 0 and n ^ 0, 
the following formula holds: 


(17) 


^ ^ /W(0)-/W(1) fi^){n) 

(27rm)^+^ (27rzn)^+^ 


Proof. The case L = 0 follows from integration by parts and the general case follows 
from iterating this result. □ 


Theorem 2. Let F he a smooth function of two real variables and assume that F 
vanishes on the boundary of R 2 , i.e., 


F{t, l-t) = F{1, t) = F{t, 1) = 0, for all t e [0,1], 

Define 

0, if {ti,t2) € Ri, 

^(^1:^2); */(^1) ^2) G i?2- 
Then the Fourier series expansion of'i’p is absolutely convergent. 


'^F{ti,t2) = 


Note: Later we will simply set </? = ik+ 1. 

Proof. We will prove this result by deriving a suitable decay estimate on | (r 1, r2) |, 

where 

^F{ri,r2)=j^ F{F,t2)e-^^"^*-dt2^dh. 


The case ri = 0 or r 2 = 0. Due to the symmetry on the assumptions on F, it 
suffices to consider only one of these cases. Let us assume that r 2 = 0. Define 

Jlifi) = [ F{ti,t2)dt2 

Ji-ti 


so that 


(18) 


^'f(ti, 0) = ii(r-i). 


Clearly we have Ti(0) = Ji(l) = 0. Setting L = 1 and / = Ji in Lemma |21 we see 
that for ri 7^ 0 

(19) |^F(ri,0)| < (^1 J((0) - J((l)| + I Jri) = O . 

With a similar estimate for |'iE'i^’(0,r2)|, we have 

X] (I’I'^’(dO)I + |4'F(0,r)|^ < 00 . 

r/0 
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The case ri ^ 0 and r 2 0- We will derive a general formula for d'F(ri,r2). To 
do this, we momentarily forget that F vanishes on the boundary of R 2 , and for 
t G [0,1], define the following functions: 


9o{t) = 

F{t,l- t), 

ho{t) = 


9i{i) = 

{d2F){t,l- t), 

hi{t) = 

{d2F){l,t), 

52 (t) = 

{did2F){t, 1-t), 

h2{t) = 

idid2F){t,l), 

53(i) = 

{did^F)(t,l-t), 

hoit) = 

(5i9|F)(l,t). 


We start with the formula for ^F{ri,r 2 ) above. Integrating by parts in the second 
variable, we obtain 


'^Firi,r 2 ) = dhe 


— 2'Kiriti 


^2-wir2t2 


- - F{ti, t2) 


*2 = 1 


-27rir2 Jt2=i-t 

p — 2Tzir2t2 

- dt2 - — - {d2F){tiR2) 

7l-ti -2TTir2 


{2TTi)r2 


(^goin - r 2 ) - ho{ri) + #a2F(*’i, *’2)) . 


We apply the same method to *’2), but integrate by parts in the first 

variable. This results in 

1 


^F{ri,r2) = “ ^ 2 ) - h 

+ “ ^ 2 ) - hi{r2) + ^did2F{ri,r2)) 

[27n)^rir2 V / 

We repeat the first two steps in the same order, which gives us 

- ^2) - ^0 


( 20 ) 


{2Fi)^rir2 

1 


(51(^1 - *’2) - lii 
(52 (?'i -r2)-h2 


(27ri)3rir2 

+ {2FiYrlrl + ^dldlF{rur2)) 

Note that from our assumptions on F, we have go = ho = hi = 0. We will have 
two subcases: 

(1) ri = r 2 = r. In this case, we easily see from the second formula above that 


(21) 


|^F(?’,'r)| < 


|gi(0)| + II^OiasFlIi 

47^2^2 


(2) ri 7^ r2. This case is slightly more subtle. We first note that gi(l) = 
/ii(0) = 0. It is also true that gi(0) = (92T’)(0,1). To see this, note that 
(9iF)(0,1) = 0 and VT'(0,1) • (1,-1) = 0, both of which follow from the 
fact that F vanishes on the boundary of i?2. The function gi being smooth 
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otherwise, we conclude that 


l5i(^i - r2)\ < 


K(l)-5i(0)l 


dir^lri — r2p 

The estimates for g 2 and /12 are simpler in nature. We use the bounds 

|52(1)-52(0)1+ I|5^||i 


|52(ti - r2)| < 


27r|ri - r2| 


as well as 




\h2{l)-h2m + \\h'2\\l 

27r|ri| 


| 53 (?’i - rz)! < II53II1, 

|^3(ri)| < llhalli, 


and 

Putting all these together, we see that 


I^f(p1)T2)| ^(i I/ m2~^22) 

V |Tir2|(ri - r2)^ |Ti(ri - r2)|r2 rfr^J 

which is easily verified to be summable over all admissible values of ri and 
r2. We will return to this shortly for a more explicit estimate. 

□ 


Explicit numerical estimates. In this subsection we will work with the specific 
function tp in (d for which 

F(ti,t2) = -40(l-ti)(l-t2) (1- (2-ti -t2)®) . 


Estimating the value o/Caxiai- Since F is symmetric, we have (p{r, 0) = (,3(0, r). We 
use the formula m to evaluate (,3(ri,0). It is a simple calculation to show that 

Ji(ti) = -15(1 - ti) + ^(1 - - 20(1 - hf + ^(1 - hf. 

Using this expression, we find that | J((0) — J((l)| = 15 and 

Hence by d, we obtain the estimate 


If we define 


|(^(r,0)| 


1^(0, t)| < 


15 + 8715 1 

471^ 


(22) axiai(iV)= E (l+(u0)| + 7(0,r)|), 

\r\<N 


then it follows that 


0 < Uaxial — C'axial(iV) < 


15 + 8715 1 
TT^ N 
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-e—B—B—B—B • 




• H—B—B—B—B- 


Figure 2. 

To estimate CaxiaK-^), we still need the actual expression for c^(r,0), which is 
given in ii^ . Taking N = 50000, numerical computation shows that Caxied{N) = 
2.90278 ..hence it follows that 

2.90278 < Caxiai < 2.90289. 

Estimating the value o/Cmain- We shall estimate the diagonal terms first. We have 

5i(t) =-240t(l-t) 

from which we obtain 

l5i(0)| = 40, 

and 

\did2F{tut2)\ = |l200(l-ti)(l-t2)(2-ti-t2)^ + 280(2-t2)®-40|, 

< 1200(1 - ti)(l - t2)(2 -h- t2)‘^ + 280(2 -ti- t2f + 40 

from which we obtain 

ll'f'ai92F||Li([o,i)2) < 80. 

Hence by we obtain the estimate 


30 1 










































A NEW UPPER BOUND FOR FINITE ADDITIVE BASES 


17 


It follows that 


i: 


|r|=A+l 

We next estimate (^(ri,r2) in the case when ri ^ r2. We begin by noting that 

120 


51(71 - 72) = 


7r2(ri - r2)2 ’ 


7l ^ 72. 


We have 


from which we obtain 

52(71 - 72) 


52(0 = 240(1 + 5t(l-t)) 
hi{t) = -40 + 280(1 -t)® 

53(0 = 240(-12- 15t + 20t^) 

= -1680(1 -t)® 

14400(2 - ti - t2)^(5 + tl + tl- 5ti - ht2 + 3tit2) 


600 


7r2(ri - r2Y ’ 
280 


7l ^ 72. 


1*2(71)1 < 

7r|ri| 

153(71-72)1 < 3080 

1*3(71)1 < 280 
ll^ofa|Flli = 2800. 

Putting these together, we obtain the estimate 

105 1 


7l ^ 0. 


(23) 


|+(7i,72)| < 


420 1 

TT^ |7ir2|(7i — r2)^ rlrl 


The following is a simple lemma: 

Lemma 3. 


(24) 


(25) 


E 


E E 

-AT+l 

E 


1 


j. 2^2 

R=Ar+l maxCT-i I ,|r 2 l) = R 1 2 

min(|ri | ,|T.2l)7iO 


< 


47r^ 1 

'Yn 


^ 4 

|7i72|(7i — r2)^ I 3 

R=N+1 max(|ri|,|r2l) = R I ^ i \ 

min(|T.i|,|r2l)#0 

ri^r2 


1 

N 


2 

max(|ri|,|r2l) = R ^ 2 


1 _ 8 1 ^ 
+2 J^2 J.2 — 

r—1 


Proof. The first inequality simply follows from 

R-l 

E 

ri|,|r 

min(|ri|,|r2l)7^0 

For the second inequality, we first use the symmetries to write 

R 


47r^ 1 

~~Y' 


E 


_ 4 

~ ~R 2^ 


1 


|7i72|(7i — 72)^ R ^ r{R — rY 

max(|ri|,|T.2l) = R I J- ^ I V J- 1 V / 

min(|ri|,|r2l)5^0 

ri#r2 


:E: 


R ^ r(R + rY 

r—1 


1 













18 


C. SINAN GUNTURK AND MELVYN B. NATHANSON 


Using the identity 


1 


1 


1 


r[R — r)2 Rr{R — r) R{R — r)^ 

and Cauchy-Schwarz inequality we have 


. R-l 

-y 

R 


1 






1 


R-l 


E l 1 “IVl 1 

(R ^ 


471^ 1 


R r(R- r)2 R2 \ ^ r(R - r) Z^ (R - 

r—1 ^ ' \r—l ^ ' r—1 ^ 

For the remaining terms, we use the trivial estimate 


r)2 


3 i?2' 


:E' 


1 


1 4 

R^ riR + r)2 “ ^ ^ 

r—1 ^ ' 


Hence 


_^_ <4f^ 

i\ M h R l’'i’’2|(ri - r 2)2 \ 3 

max(|rx|,|r2l) = -R I I v / V 

min(|ri |,|r 2 |)?i0 


E 


i?2 


and the result follows. 


If we define 
(26) 


N 


a 


.(«) = E E 


\‘f{ri,r2)\, 


R—1 max(|ri|,|r2l) = R 

min(|ri|,|r2l)7^0 


□ 


then we have 


(27) 


0 E lEniain ~ C'niain(.^) < 



420\ 1 40 

TV ^ iV' 


For N = 4000, numerical computation using the formulas and m reveals 
that C'niain(.^) = 4.75145 ...; hence with the above error estimate, we have 


(28) 


4.75145 < Cmain < 4.76146. 


Explicit expressions for ‘f{ri,r 2 )- The following formulas have been computed 
using Mathematica, though it is also possible to compute them easily using the 
iterative procedure based on integration by parts which was outlined in this section 
earlier. 


(29) (^(r,0) 


47r2r2 y 7r2r2 


45 135 \ 

TT^r® J 

. 60 / 63 315 945 \ 

167r®r®y 


N 10 21 315 945 \ 

(30) + 

. 55 / 126 


630 945 'N 

Il7r4^4 llTr^r^J 
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(31) (^(r,s) = 

1575 


525 


35 


1575 


47r®r'®(r —s) 47r®r^(r —s) 2 (r — s) 47r® (r—s) s® 


+ 


225 


525 


225 


+ 


225 


2 TT® r (r — s)^ s® 4 tt® (r — s)^ 2 tt® (r — s)^ 2 tt® r® (r — s)^ s® 


75 


35 


225 


75 


2 7r®r(r-s) s® 2 7r4(r-s) 2-K^r'^{r-s) 2 7r®r2(r-s) s2 

225 75 5 


+ 


+ 


2 7r®r5(r-s)^s 2 tt® r® (r - s)^ s r {r - sf s 


4-^ 

V 4 TT^ 


1575 


525 


105 


1575 


47r9r7(r-s)^ 2'k'^ {r - sf 2TT^r^{r-sf 47r9(r-s)^s7 


+ 


225 


525 


225 


+ 


225 


2 7r9r(r-s)^s® 2 7r7(r-s)^s5 2 tt® r2 (r - s)^ s® 2 tt^ r® (r - s)^ 


75 


105 


225 


75 


+ 


7r7r(r-s) 2^^ {r - s) s® 2 7r9r4(r-s) s® {r - s) s® 

225 75 15 225 


2 7 r®r®(r —s) s2 tj- 7 j-3 _ g'j g 2 jj-5 j. _ g'j g 2 2 7 r®r®(r —s) s 

75 15 \ 


r7 /r’4 


(r - s)^ 


S TT® r2 


(r - s)^ 
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